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Abstract 

We report an empirical study of the Ibovespa index of the Sao Paulo Stock Exchange 
in which we detect the existence of long-range correlations. To analyze our data we 
introduce a rescaled variant of the usual Detrended Fluctuation Analysis that allows 
us to obtain the Hurst exponent through a one-parameter fitting. We also compute 
a time-dependent Hurst exponent H{t) using three-year moving time windows. In 
particular, we find that before the launch of the Collor Plan in 1990 the curve H{t) 
remains, in general, well above 1/2, while afterwards it stays close to 1/2. We thus 
argue that the structural reforms set off by the Collor Plan has lead to a more 
efficient stock market in Brazil. We also suggest that the time dependence of the 
Ibovespa Hurst exponent could be described in terms of a multifractional Brownian 
motion. 

Key words: Long memory processes, Dentrended fluctuation analysis. Hurst 
exponent, Econophysics, Multifractional Brownian motion 
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1 Introduction 



In recent years, it has been realized that many problems from Economics can 
be studied with the standard "toolkit" of statistical physics [1,2]. For instance, 
asset prices in financial markets are commonly described in terms of a geo- 
metric Brownian motion, an assumption that is at the heart of the so-called 
Efficient Market Hypothesis (EMH). Thus, in the EMH scenario the returns 
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on a given stock follow an uncorrelated Gaussian process (white noise). In 
more common parlance, the EMH says that all information available about a 
given financial asset is already reflected on its current price, so that knowing 
the asset past history does not in any way help us to predict future prices. 
Although the EMH stands as a cornerstone of modern Finance [3] , deviations 
from efficiency have been recently observed in many different financial mar- 
kets [1,2]. In such 'inefficient' markets the empirical data violate cither the 
independence or the Gaussian assumptions of the EMH. The former instance 
of 'inefficiency' is of particular interest to us here because it imphes the ex- 
istence of long-range correlations that are not accounted for in the standard 
EMH model. 

A time series can be tested for correlation in many different ways [4]. One 
general methodology consists in estimating how a certain fluctuation measure, 
to be denoted here generically by F, scales with the size r of the time window 
considered. Specific methods, such as the Hurst rcscaled range analysis [5] or 
the Detrendend Fluctuation Analysis [6,7], differ basically on the choice of the 
fluctuation measure; see below for more details about these two methods. If 
the time series is uncorrelated one expects that F ~ r^/^, as is the case for the 
standard Brownian motion. On the other hand, if F ~ with H ^ 1/2 one 
then says that the time series has long-term memory, with if > 1/2 [H < 1/2) 
meaning persistence (antipersistence). The exponent H is generally referred 
to as the Hurst exponent. 

In the present paper we analyze the behavior of the Sao Paulo Stock Exchange 
Ibovespa index, which is the main index of the Brazilian stock market. We 
have carried out a Dctrended Fluctuation Analysis (DFA) of the Ibovespa 
covering over 30 years of data, since its inception in 1968 until the year 2001. 
Here we introduce however a rescaled variant of the original DFA [6] that 
has the advantage of allowing us to determine the Hurst exponent H with 
a one-parameter fitting. For the complete Ibovespa time series we find H — 
0.6, indicating that the Ibovespa index exhibits persistence. We have also 
computed a 'local' Hurst exponent H{t) in moving three-year time windows 
and found that H(t) varies considerably over time. In particular, we find that 
the curve H{t) undergoes a distinct change in character around the year 1990. 
More specifically, we observe that before 1990 the exponent H{t) is always 
greater than 1/2. Then around the year 1990 H{t) drops quite rapidly towards 
1/2 and stays (within some fluctuation) around this value afterwards. We 
identify this drastic change in H{t) as a consequence of the economic plan 
adopted in March 1990, the so-called CoUor Plan, that marked the beginning 
of structural reforms in the Brazilian economy. This interpretation is confirmed 
by separate analyses of the Ibovespa returns prior and after the Collor Plan, 
where in the former case we find H = 0.63 while for the latter period we have 
H pa 0.5. We thus argue that the process of opening and modernization of the 
economy started by the Collor Plan has led to a more efficient stock market 
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in Brazil, in the sense that the Hurst exponent H for this period is close to 
the value predicted by the EMH. 

The fact that its Hurst exponent changes over time indicates that the Ibovespa 
follows a multifractal (rather than monofractal) process. We thus suggest that 
the Ibovespa basic dynamics could perhaps be captured by the so-called mul- 
tifractional Brownian motion (MFBM) [8], which is a generalized fractional 
Brownian motion with a time-dependent Hurst exponent H{t). One interest- 
ing property of the MFBM is that it is a Gaussian additive process, while the 
multifractal processes often considered in physics [9,10] and economics [11] are 
in general multiplicative. 

The paper is organized as follows. In Sec. 2 we describe our data. In Sec. 3 we 
present some background material about fractional Brownian motion and give 
a brief description of the Detrended Fluctuation Analysis. In particular, we 
discuss a rescaled version of the DFA that allows us to determine the Hurst 
exponent H through a one-parameter fitting. The results of this rescaled DFA 
for the Ibovespa time series are shown in Sec. 4. In Sec. 5 we offer a succinct 
discussion of the multifractional Brownian motion and its possible relevance 
to financial time series. In Sec. 6 we review our main results and conclusions. 
For completeness, we present in Appendix A the results of the Hurst analysis 
for our Ibovespa time series. It is shown there that in this case the Hurst 
method is shghtly less reliable than the DFA in detecting long- correlation, 
thus justifying our choice of restricting ourselves to the DFA in the main 
portion of the present paper. 



2 The Data 



The Ibovespa index represents the present value of a self-financing hypothetical 
portfolio made up of the most traded stocks on the Sao Paulo Stock Exchange. 
In this paper we have analyzed the daily closing values of the Ibovespa since its 
inception (January 02, 1968) until recent times (May 31, 2001), amounting to 
8209 trading days. The nominal figures of the daily Ibovespa were deflated by 
Brazil's general price index (IGP-DI, in its Portuguese acronym) and converted 
to Brazilian Reals (BRL) in values of August 1994. The corresponding time 
series for the deflated Ibovespa is shown in Fig. 1, while the IGP-DI price 
index is plotted in Fig. 2 in monthly rates. 

In the period comprised in our study many economic events took place in 
Brazil and abroad that had a direct impact on the Brazilian stock market. For 
the benefit of the reader less familiar with Brazil's recent economic history, it 
was thought desirable to recall here some of the main events whose effects can 
be clearly distinguished in the data shown in Figs. 1 and 2. For instance, the 
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first large peak in the Ibovespa in the early 1970's corresponds to the so-called 
'Brazilian economic miracle' (during the military regime that took power in 
1964 and ended in 1985) when Brazil's GDP grew over 10% a year. After the 
first oil crises in 1973 the Brazilian economy entered a long period of slower 
growth and steadily increasing inflation, as reflected in the consistent decline 
of the Ibovespa from 1973 to 1984 (Fig. 1) and a corresponding surge of the 
inflation rate (Fig. 2). From 1986 to 1991 six largely imsuccessful economic 
plans were adopted by the succeeding governments in attempts to control 
inflation. The effects of these economic plans are very clearly distinguished in 
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Fig. 1. Daily closing values of the deflated Ibovespa index in the period January 
1968-May 2001. 
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Fig. 2. Brazilian monthly inflation rate in the period January 1968-May 2001. The 
arrows indicate the time of launch of the following economic plans: Cruzado Plan 
(A), Collor Plan (B), and the Real Plan (C); see text. 
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Fig. 3. Daily returns of the Ibovespa in the period January 02, 1968 {t = 0) through 
May 30, 2001 it = 8207). 

Fig. 2 as they correspond to sudden drops in the inflation rate. For example, 
the first of such plans, the so-called Cruzado Plan, was launched in February 
1986 and led not only to a sharp decline in the inflation rate (arrow A in 
Fig. 2) but also to a strong peak in the Ibovespa index (see Fig. 1). Shortly 
afterwards, however, inflation picks up again and the Ibovespa recedes to pre- 
Cruzado levels. Although it also failed to rein in inflation, a plan of more 
lasting impact on the Brazilian economy was the CoUor Plan (arrow B in 
Fig. 2), as we shall see later. It should also be recalled that the Real Plan, 
launched in July 1994 (arrow C in Fig. 2), has so far succeeded in stabilizing 
the economy and keeping inflation at bay. 

In what follows we shall carry out a detailed analysis of the Ibovespa daily 
returns, which we deflne in the usual way. Let yt denote the closing value of the 
Ibovespa at time t, where t is measured in trading days and t = corresponds 
to the date January 02, 1968, when the Ibovespa was created. The return at 
time t is given by the logarithmic difference between consecutive daily values: 



In Fig. 3 we plot the daily returns for the Ibovespa index. In what follows 
we wish to investigate the existence of long-term dependence in the Ibovespa 
daily returns. Before going into this analysis, however, we shall first present 
some background material about correlated time series and describe briefiy 
the method we will use to detect correlation in our data. 
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3 Correlated Time Series 



One of the simplest and most important stochastic models that display long- 
term dependence is the fractional Brownian motion (FBM), whose main prop- 
erties we will briefly review below. Following this, wc shall then describe the 
Detrended Fluctuation Analysis (DFA) which, as already noted, is the main 
method we will use to estimate the Hurst exponent H for our empirical data. 
(Several other estimators for H have been discussed in the literature; see, e.g., 
Ref. [4] for a comparison among some of them.) 

3.1 Fractional Brownian motion 

We recall here that the fractional Brownian motion {BH{t),t > 0} is a Gaus- 
sian process with zero mean and stationary increments whose variance and 
covariance are given by 

E[5|(t)] = (2) 
E[BH{s)BH{t)] = (5^^ + - \t - s^) , (3) 

where 0<i?<l,(7>0, and E[] denotes expected value. The process 
Bffit) is statistically self-similar (or more exactly self-affine) in the sense of 
finite-dimensional distributions: 

BH{at) = a'^Bnit), (4) 
for all a > 0, where = means equality in distribution. 

The parameter H is called the self-similarity exponent or the Hurst exponent. 

For if = 1/2 the process B^it) corresponds to the standard Brownian motion, 
in which case the increments Xt = Bnit + l) — Bnit) are statistically indepen- 
dent and represent the usual white noise. On the other hand, for H ^ 1/2 the 
increments Xt, known as fractional white noise, display long-range correlation 
in the sense that 

E[Xt+hXt] ~ a^2H{2H - l)h^"-'^ for h ^ 00, (5) 

as one can easily verify from (2) and (3). Thus, ifl/2<if<l the increments 
of the FBM are positively correlated and we say that the process Bnit) ex- 
hibits persistence. Likewise, for < H < 1/2 the increments are negatively 
correlated and the FBM is said to show antipersistence. 
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3.2 Detrended fluctuation analysis 



The Detrended Fluctuation Analysis (DFA), which is a modification of the 

usual variance analysis, was proposed independently albeit with different names 
in Refs. [6] and [7] . (In Ref. [7] the analogous of the fluctuation function F{t) 
defined below was termed Straight Line Roughness.) The DFA has the ad- 
vantage over the standard variance analysis of being able to detect long-term 
dependence in nonstationary time series [6]. The idea of the method is to 
subtract possible deterministic trends from the original time series and then 
analyze the fiuctuation of the detrended data, as explained below. 

To implement the DFA, we first integrate the original time series {rt}t=i,...,T 
to obtain the cumulative time series X{t): 

X(t)^J2(n'-r), t^l,...,T, (6) 
t'=i 



where 

1 ^ 



r^Tf^Y^rt'- (7) 
^ t'=i 



Next we break up {X{t)} into N non-overlapping time intervals, /„, of equal 

size r, where n = 0, 1,...,A^ — 1 and corresponds to the integer part of T /t. 
We then introduce the local trend function defined by 

n(0=a„ + M for teln, (8) 



where the coefficients and hn represent the least-square linear fit of X{t) in 
the interval /„. We now compute the fiuctuation function F{t) defined as the 
root mean deviation of X{t) with respect to the trend function Yr{t): 



\^Y.\m-yr{t)]\ (9) 



F(r) = 



In the original DFA [6] the Hurst exponent H is obtained directly from the 
scaling behavior of the above fiuctuation function: F{t) ~ . For reasons 
that will become apparent shortly, we find it convenient to introduce a rescaled 
(dimensionless) fiuctuation function Fs{t), 

Fs{r) ^ (10) 
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where S is the data standard deviation 



(11) 



The rescaled fluctuation function Fs{t) obviously has the same functional 
form as the original F{t) and so we write 

Fs{t) = Cr^, (12) 



where C is a constant independent of the time lag r. From here on, we shall 
drop the subscript from Fg{r), with the understanding that we will be working 
solely with the rescaled fluctuation function. 

In a double-logarithmic plot the relationship (12) yields a straight line whose 
slope is precisely the exponent H, and so the empirical value for H can be 
easily obtained from a linear regression (in log-log scale) of the corresponding 
data for F{t). One practical problem with this method, however, is that one 
needs to choose an appropriate interval within which to perform the linear fit. 
For instance, the first few points at the low end of the graph of F(r) should be 
disregarded because in this region the detrending procedure removes too much 
of the fluctuation — this effect accounts for the usual 'bending down' of F{t) 
for small r. For large r, on the other hand, there are few boxes /„ for a proper 
averaging to be made and hence the values of F{t) are not statistically rehable 
in this case. Furthermore, it has been found [12] that the values obtained for 
H using this procedure are somewhat dependent on the choice of the fltting 
interval. In order to avoid some of these difficulties we shall use here a slightly 
different procedure to estimate the exponent H from the graph of F{t). 

To do this, we shall rely on the fact that for the fractional Brownian mo- 
tion, i.e., for X(t) = Bnit), the (rescaled) fluctuation function F{t) can be 
computed exactly [4]: 

Fh{t) = Cj^T^, (13) 



where 

2 12 

Ch = + . (14) 

In (13) we have added a subscript H to the fluctuation function F{r) to de- 
note explicitly that it refers to Bfjit). Now, to estimate H for our time series 
we shall simply adjust the parameter H so as to obtain the best agreement 
between the theoretical curve predicted by Fh{t) and our empirical data for 



8 



F(t). In this way, we can estimate H with a one-parameter fitting! In princi- 
ple, one could carry out a formal nonlinear regression using (13) and (14) to 
determine H, but in practice it suffices to vary H incrementally and decide by 
visual inspection when the theoretical curve best matches the empirical data. 
We shall now apply this methodology to estimate the Hurst exponent for the 
Ibovespa returns. 



4 Fluctuation analysis of the Ibovespa 

In this section we discuss the results of the DFA applied to our Ibovespa time 
series, which we recall spans the period from January 1968 through May 2001. 
In Fig. 4 we plot in double-logarithmic scale the corresponding fluctuation 
function F{t) against the window size r. Using the procedure outlined at the 
end of the previous section, we obtain the following estimate for the Hurst 
exponent: H = 0.60 ± 0.01. In Fig. 4 we have also plotted (straight line) the 
corresponding theoretical curve F/y(r) given by (13) and (14), which as one can 
see is in excellent agreement with the empirical data for intermediate values 
of r. Since H > 1/2 we conclude that the Ibovespa returns show persistence. 
We also note that for r > 130 (indicated by the arrow in Fig. 4) the empirical 
data deviate somewhat from the initial scaling behavior and appear to cross 
over to a regime with a slope closer to 1/2. This indicates that the Ibovespa 
tends to loose its 'memory' after a period of about 6 months. It is perhaps 
worth pointing out that similar behavior was seen, for example, in several 
equity indices of the London Stock Exchange [13]. There, the break from the 
scaling regime typically occurred around 160 trading days [13], which is not 
very different from what we found for the Ibovespa, especially considering that 
the location of the crossover point is difficult to determine precisely. 

In order to check the validity of our analysis above, we have also computed the 
Hurst exponent for a shuffled version of our data in which we have randomly 
mixed the time series of the Ibovespa returns. Because the process of shuffling 
tends to destroy any previously existing correlation, we would now expect a 
Hurst exponent equal (or very close) to 1/2. This is indeed the shown 
in Fig. 5, where we plot F{t) for the shuffled data. In this figure one sees 
that the data points are extremely well described by the theoretical curve 
(straight line) given by Fi{{t) with H = 0.5. This thus confirms the fact that 
the long-term dependence seen in the original time series was not an spurious 
effect. 

We now wish to investigate whether the Hurst exponent for the Ibovespa varies 
in time, which would indicate the existence of nonstationary fluctuations in the 
Brazilian stock market. In fact, a simple visual inspection of the distribution of 
the Ibovespa returns shown in Fig. 3 already revels evidences of nonstationary 
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Fig. 4. Fluctuation function F{t) as a function of window size r for the returns of 
the Ibovcspa index in the period 1968-2001. The straight Une gives the theoretical 
curve Fh{t) for H = 0.6; see text. 
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Fig. 5. Fluctuation function F{t) for the shuffled Ibovespa returns. The straight 
line represents the theoretical curve Fh{t) for H = 0.5 

behavior, which we would like to quantify better. Furthermore, in the recent 
past not only the Brazilian economy was plagued by runaway inflation and 
but it also endured several 'shock-therapy' economic plans, as we have already 
mentioned. It is thus natural to ask how such stressful events may have affected 
the dynamics of the Ibovespa index, particularly in regard to its degree of 
correlation. We now turn our attention to this point. 

In order to estimate a time- varying Hurst exponent using a fluctuation analysis 
(or any method for that matter) one is presented with a difficult challenge, 
for one has to try to satisfy competing requirements. Indeed, to determine a 
'local' exponent if at a particular time t one has to consider a time interval 
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Fig. 6. The Hurst exponent H for the Ibovespa as a function of time. Here H was 
computed in three-year periods and the variable t represents the origin of each such 
interval. 

around t that is considerably smaller than the total time spanned by the data 
but still sufficiently large to contain enough points for a meaningful statistics. 
In the case of the Ibovespa we have found that a three- year period (736 trading 
days) is an acceptable compromise between these two opposing demands. (For 
smaller periods the fluctuations in the values of H are too large to be trusted.) 
We have accordingly applied the DFA to the Ibovespa returns in moving three- 
year intervals. More specifically, our methodology is as follows: starting at the 
beginning of our time series we compute the exponent H considering only the 
data points within three years from the initial point, then we advance our 
three-year time window by 20 days and compute H for this new period, and 
so on. Proceeding in this way we obtain the curve for H{t) shown in Fig. 6, 
where t denotes the origin of each three-year time interval. 

The most striking feature seen in Fig. 6 is perhaps the fact that the Ibovespa 
shows persistence {H > 1/2) all the way up to the early 1990's, after which it 
'switches' to a regime with alternating persistent and antipersistent behavior 
but where H remains somewhat close to 1/2. Notice in particular that during 
the 1970's and 1980's the curve H{t) stays well above 1/2, the only exception to 
this trend occurring around the year 1986 when H dips momentarily towards 
1/2 — we associate this effect with the launch of the Cruzado Plan in February 
1986; see Sec. 2. Shortly after the Cruzado Plan, however, infiation picks up 
again (Fig. 2), the Ibovespa declines (Fig. 1), and the exponent H returns 
to prc-Cruzado levels (Fig. 6). Then in the early 1990's, at the time of the 
CoUor Plan (see below), we observe a dramatic decline in the curve H{t) 
towards 1/2, with H remaining (within some fiuctuation) around this value 
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Fig. 7. Plot of F(t) for the Ibovespa returns before and after the CoUor Plan. The 
upper (lower) straight line corresponds to the theoretical curve Fh{t) with H = 0.63 



afterwards. Note also that around the time of the launch of the Real Plan in 
July 1994 the curve H{t) reaches its lowest values. We have thus seen that the 
Ibovespa Hurst exponent declines following the adoption of a major economic 
plan (such as the Cruzado, Collor and Real Plans), thus confirming the fact 
[15] that a firm intervention on the market tends to reduce H momentarily. 

As we already mentioned, the drop in H that occurs in the early 1990's can 
be unambiguously traced back to the launch, in March 1990, of the Collor 
Plan which marked the beginning of structural reforms in Brazil [14]. In this 
context, it is important to realize that this decline in H takes place well 
before the adoption of the Real Plan. It is indeed quite surprising that the 
Hurst exponents found in the early 1990's, in a macroeconomic environment 
dominated by very high inflation, are on about the same level as those found 
after the Real Plan when the economy had been stabilized. This fact in itself is 
a clear evidence of the profound impact that the structural reforms set off by 
the Collor Plan had on the Brazilian economy. In order to study this effect in 
more detail, we have apphed the DFA separately to the Ibovespa returns prior 
and after the Collor Plan. Our results are shown in Fig. 7. In this figure one 
sees that before the Collor Plan the Brazilian stock market shows a significant 
amount of persistence {H = 0.63), while after the Collor Plan the Ibovespa 
displays essentially no long-term dependence {H = 0.5). 

In light of the preceding results, it thus seems legitimate to conclude that 

the opening and consequent modernization [14] of the Brazilian economy that 
begun in the early 1990's has led to a more efficient stock market, in the sense 
that H = 0.5 for the Ibovespa after 1990. By the same token, we associate 
the period of higher H before 1990 with a less efficient market. In those years. 



{H = 0.5). 
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the Brazilian economy was considerably closed to foreign competition and its 
financial markets were not readily accessible to international investors. In such 
closed economic environment the Brazilian stock market was conceivably more 
prone to 'correlated fluctuations' (and perhaps also more susceptible to being 
pushed around by aggressive investors), which may explain in part a Hurst 
exponent greater than 0.5 for that period. 



5 Multifractional Brownian motion 

We have seen above that the Hurst exponent H for the Ibovespa stock index 
changes appreciably over time. In this regard, the fractional Brownian motion 
appears to be a somewhat restrictive model, unable to capture more fully 
the complex dynamics of the Brazilian stock market. Similar time-dependent 
Hurst exponents have also been observed in other financial markets [15,16,17] 
and in physical processes such as traffic flow in the Internet [18] . To circumvent 
this shortcoming of the FBM, Peltier and Levy-Vehel [8] have proposed a 
multifractional Brownian motion (MFBM) in which the scaling exponent H 
is allowed to vary in time. For a rigorous exposition of the definition and 
mathematical properties of MFBM the interested reader is referred to Ref. [8]. 
Here we shall present only the main properties of the MFBM that are of 
interest to us. Before we proceed, it should however be emphasized that our 
aim in this section is not to develop a specific model for the Ibovespa but 
rather to present a generic theoretical framework based on the MFBM in which 
the time-dependence of the Ibovespa Hurst exponent could, in principle, be 
understood. 

Let H{t) be a Holder-continuous function in the interval t G [0, 1] with Holder 
exponent > 0, such that for any i > we have < H{t) < min(l, /3) [8]. (For 
ease of notation, we shall sometimes indicate the t-dependence as a subscript 
and write H^.) The multifractional Brownian motion {WHt{t),t > 0} is the 
Gaussian process defined by 

1 'r 

WhM = J at - - {-s)f-''']dB{tl (15) 

where V{x) is the gamma function, (,x)+ equals ,x if x > and zero otherwise, 
and B{t) denotes the usual Brownian motion. One important aspect of the 
process defined above is that its increments are no longer stationary since it 
can be shown [8] that 

E[{WH,^St + h) -WH,{t)Y]'^h'^"^^^ as /i ^ 0. (16) 
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Because of its non-stationarity the MFBM is no longer a self-similar process 
either. However, it is possible to define the concept of locally asymptotically 
similarity at the point to > in the following way 



Wut^ {t) as a ^ 0. 



(17) 



[Compare this with the global self-similar property of the FBM given in (4).] 
In this sense, WHt{t) can be thought of as a process that locally at time t 
'resembles' a fractional Brownian motion Bntit)- In fact, a practical way to 
generate a mult if r actional Brownian path consists in generating for each time 
t e [0,1] a fractional Brownian motion Bh^{s), s G [0,1], and then setting 
WHt{t) = BhX^)- More details about this method together with an algorithm 
for its implementation can be found in Ref. [8]. (We remark parenthetically, 
however, that the algorithm presented there for the computation of H{t) is 
not appropriate for empirical time series, as inadvertently used in Ref. [17], 
because such algorithm requires previous knowledge of the variance of the 
FBM used to generate the MFBM, which of course is not known for empirical 
data.) 

In Fig. 8 we show an example of a path of multifractional Brownian motion 
characterized by the Hurst function 



Although somewhat arbitrary, this example was chosen to mimic the generic 
trend seen in the Ibovespa where H{t) can be viewed as having (on aver- 
age) two distinct 'plateaus,' one before and the other after 1990, with a swift 
crossover between them. For comparison purpose, we have generated a MFBM 
path with approximately the same number of data points (A^ = 2^^ = 8192) 
as in the original Ibovespa time series {N — 8209). We have then apphed 
the detrended fluctuation analysis to our MFBM path, with the result being 
shown in Fig. 9. Since a MFBM path is not characterized by a constant Hurst 
exponent, the functional dependence of the DFA fluctuation function F{t) 
will not, in general, be a power law. Nevertheless, we see from Fig. 9 that 
F(r) exhibits an approximate scaling regime, F{t) ~ with H — 0.6, up 
to about r = 160. For comparison, we also show in Fig. 9 the corresponding 
relationship (straight line) predicted by Fh[t) given in (13), which as one sees 
is in good agreement with the data in the scaling region. 

We thus interpret the scaling exponent H as a, kind of average Hurst exponent 

for our MFBM time series. Note, however, that although the region of small 
H represents only about 20% of the total time series it has nonetheless a 
significant weight on the exponent H, since its value is considerably smaller 
than what one would obtain by simply averaging the curve H{t) [in this case 



H{t) = 0.63 - 0.076 arctan(30t - 24). 



(18) 
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Fig. 8. Path of a multifractional Brownian motion (a) with the Holder function 
H{t) = 0.63 - 0.076 arctan(30i - 24) shown in (b). 

one gets ifave = 0.68]. This result is in agreement with the known fact [12] 
that regions of large fluctuations (and hence smaller H) tend to dominate 
the behavior of F{t) for small-to-intcrmcdiate values of r. A more detailed 
discussion about the detrended fluctuation analysis of MFBM paths will be 
presented elsewhere. 

The preceding discussion thus suggests that the Hurst exponent H obtained 
by applying the DFA (and similar fluctuation analyses [4]) to flnancial time 
series might perhaps capture only an average behavior over the time period 
spanned the data. This seems to be particularly true for the Ibovespa index, 
where we find that the DFA Hurst exponent obtained using the complete time 
series (spanning over 30 years of data) is H = 0.6, whereas during this period 
the local exponent H{t) varies considerably within the range (0.76,0.42); see 
Fig. 6. One is thus lead to conclude that when dealing with long financial time 
series, a more localized (in time) analysis is necessary to determine a possible 
dependence of H on time. In particular, note that a global Hurst exponent 
equal to 1/2 docs not necessarily imply the absence of correlation: there could 
be positive and negative correlations at different periods of time that average 
out so as to yield an effective if = 0.5. 



6 Conclusions 

We have carried out a Detrended Fluctuation Analysis of the daily returns 
of the Sao Paulo Stock Exchange Index (Ibovespa) , spanning over 30 years of 
data from January 1968 up to May 2001. For this time series we have obtained 
a Hurst exponent greater than 1/2, indicating that the Ibovespa has long-term 
dependence (persistence). This memory effect seems to last for up to about 
6 months (130 trading days), after which time the fluctuation function F{t) 
deviates from the initial scaling behavior and crosses over to a regime with a 
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Fig. 9. Detrended fluctuation analysis of the multifractional Brownian motion shown 
in Fig. 8(b). 



slope closer to 1/2. Similar behavior has been observed in other stock indices 
[13], and it thus appears that stock markets tend to loose their 'memory' 
typically in about half a year. 

We have also performed a more localized (in time) analysis of the Ibovespa 
returns, where we have calculated the local Hurst exponent H{t) in three-year 
moving time windows. Here the most striking feature is the fact that after 
1990 the Hurst exponent stays close to 1/2 (within some fluctuation), whereas 
in the preceding decades it was, in general, considerably larger than 1/2. We 
have thus argued that the liberalizing measures (e.g., substantial reduction 
of import tariffs, cutting down the 'red tape' for foreign investments, and 
privatisation of state companies [14]) brought about by the CoUor Plan in 
the early 1990's and continued with the Real Plan in 1994 resulted in a more 
efficient stock market in the sense that H remained close to 1/2 after 1990. 

The fact that its Hurst exponent changes over time time indicates that the 
Ibovespa time series cannot be satisfactorily modeled in terms of a simple 
fractional Brownian motion. In other words, the Ibovespa follows a multifractal 
process. We have thus suggested that the multifractional Brownian motion, in 
which the scaling exponent H is allowed to vary in time, could perhaps be a 
more appropriate description for the Brazilian stock market dynamics. It thus 
remains an interesting question to investigate to what extent the MFBM is a 
valid model for other financial time series. 
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A Appendix: Hurst rescaled range analysis 



The Hurst rescaled range {R/ S) analysis is a technique devised by the hy- 
drologist Henry Hurst in 1951 [5] to test for the presence of correlations in 
empirical time series. A detailed description of the Hurst analysis (including 
its historical origin) can now be found in textbooks [9,10], so that here we 
shall give only a brief summary of the method and then proceed to apply it 
to the Ibovespa returns. 

In the Hurst analysis, we start by dividing our time series {rt}t=i,...,T in N non- 
overlapping time intervals /„ of equal size r. For each interval we compute 
the range 



Rn = max 

l<fe<r 



i=l 



mm 

l<fe<r 



i=l 



(A.l) 



where 



t=l 



(A.2) 



and the data standard deviation 



T 



(A.3) 



t=i 



The rescaled range (R/S)-,- is then defined as the average of the ratio Rn/ 
over all intervals 
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Fig. A.l. Rescaled range R/S versus the time lag r for the Ibovespa returns in 
the period 1968- 2001 and for the corresponding shuffled data. Also shown for 
comparison is a straight line with slope equal to 1/2. 

The Hurst exponent H is obtained from the scaling behavior of {R/S)t: 



A plot of the rescaled range R/S as a function of r for the Ibovespa returns 
over the entire period (1968-2001) analyzed in this paper is shown in the upper 
curve of Fig. A.l. The data in this case show a scaling regime that goes from 
T = 10 up to about r = 160. A linear regression in this region yields the value 
H — 0.65, which is somewhat larger than the exponent {H — 0.6) obtained 
via the DFA for the same time series. In order to understand the origin of such 
discrepancy, we have also applied the Hurst method to the shuffled Ibovespa 
data (lower graph in Fig. A.l) and found H = 0.57 in this case. We thus see 
that the Hurst analysis of the shuffled data predicts an unexpected residual 
correlation. (Recall that applying the DFA to the shuffled data yields H — 0.5, 
as expected.) This 'excess of correlation' found in the R/S analysis of the 
shuffled Ibovespa is perhaps not entirely surprising, given that it is known 
that the Hurst method tends to overestimate the Hurst exponent for time 
series of small sizes [9]. This effect may perhaps explain why the exponent H 
obtained via the Hurst method is usually larger than that found via the DFA 
[19]. We thus conclude that the DFA appears to be a more reliable method 
for detecting true correlation in financial time series. It is worth therefore 
emphasizing that the numerical values obtained for the exponent H via the 
Hurst method should be interpreted with some caution, specially when the 
time series is not sufficiently large, as is commonly the case in practice. 



(R/S)^ ~ r^. 



(A.5) 
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